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ABSTRACT 

The observational evidence that Super-Massive Black Holes (M, ~ 10^“are already in 
place less than 1 Gyr after the Big Bang poses stringent time constraints on the growth effi¬ 
ciency of their seeds. Among proposed possibilities, the formation of massive (~ 10^“® M©) 
seeds and/or the occurrence of super-Eddington (M > MecIci) accretion episodes may con¬ 
tribute to the solution of this problem. In this work, using a set of astrophysically-motivated 
initial conditions, we analytically and numerically investigate the accretion flow onto high- 
redshift (z~ 10) black holes to understand the physical requirements favoring rapid and effi¬ 
cient growth. Our model identifies a “feeding-dominated" accretion regime and a “feedback- 
limited" one, the latter being characterized by intermittent (duty cycles D < 0.5) and inef¬ 
ficient growth, with recurring outflow episodes. We And that low-mass seeds (< 10^“^M©) 
evolve in the feedback-limited regime, while more massive seeds (> 10^“®M©) grow very 
rapidly as they are found in the feeding-dominated regime. In addition to the standard ac¬ 
cretion model with a Axed matter-energy conversion factor (e = 0.1), we have also explored 
slim disk models, appropriate for super-Eddington accretion, where radiation is trapped in the 
disk and the radiative efficiency is reduced (e < 0.04), which may ensure a continuous growth 

with M » MecIcI (up to ~ 3Q0MEdd in our simulations). Under these conditions, outflows play 
a negligible role and a black hole can accrete 80% - 100% of the gas mass of the host halo 
(~ 10^M©) in ~ lOMyr, while in feedback-limited systems we predict that black holes can 
accrete only up to ~ 15% of the available mass. 

Key words: accretion - black hole physics - quasars: supermassive black holes - cosmology: 
theory - dark ages, reionization, first stars - cosmology: early Universe 


1 INTRODUCTION 

Recent observations (Mortlock et al. 2011; Wu et al. 2015) have de¬ 
tected the presence of optically bright quasars at redshifts as high 
as z ~ 7. These high-energy sources are powered by accretion onto 
Super-Massive Black Holes (SMBHs), suggesting the presence of 
compact objects with mass M, ~ lO^^'^’M© less than IGyr after 
the Big Bang (Fan et al. 2006). This evidence contrasts with the 
standard theory of black hole growth, which requires a longer time 
to build up such massive objects (see Haiman 2013 for a recent re¬ 
view), due to: (i) the low mass of some of the proposed seeds, bom 
out of first-generation (Pop III) stars, with masses M, < 10^ M© 
(Madau & Rees 2001; Bromm & Loeb 2003; Petri et al. 2012), 
and (ii) the maximum growth rate allowed for a radiatively efficient 
and spherical inflow, the Eddington rate, which provides a lower 
limit for the time scale of the process (Jeon et al. 2012). Generally 
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speaking, the luminosity L emitted due to a gas inflow with an ac¬ 
cretion rate M = dMjdt is L = ec^M, where e is the matter-radiation 
conversion factor and c is the speed of light. In the standard sce¬ 
nario, a black hole grows in mass exponentially, with a time scale 
given by the Salpeter time ts ~ 0.045eo.l Gyr, where eo.i is normal¬ 
ized to the standard value of 10%: starting from a low-mass seed 
(M, ~ 100Mq), this process would require a constant accretion at 
the Eddington rate up to z ~ 7 to build up a ~ 10® M© SMBH. 

Several ways to overcome these limitations have been pro¬ 
posed in the literature. The possibility of giving a jump start to the 
growth process through more massive seeds has been investigated 
thoroughly (see e.g. Volonteri 2010) with a variety of mechanisms: 
(i) the direct collapse of self-gravitating pre-galactic disks at high- 
redshifts (Lodato & Natarajan 2006; Begelman et al. 2006; Lodato 
& Pringle 2007), (ii) the formation of a very massive star from mn- 
away stellar mergers in a dense cluster (Devecchi & Volonteri 2009; 
Davies et al. 2011) and (iii) the collapse of a primordial atomic- 
cooling halo, exposed to a Lyman-Werner flux of sufficient inten- 


© 0000 The Authors 


2 F. Pacucci et al. 


sity, into a Direct Collapse Black Hole (DCBH), through a general 
relativistic instability (Shang et al. 2010; Johnson et al. 2012; Fer¬ 
rara et al. 2014). 

An alternative scenario assumes that accretion rates are not 
capped by the Eddington limit (Volonteri & Rees 2005). Recent 
works (Alexander & Natarajan 2014; Madau et al. 2014; Volonteri 
et al. 2014) have proposed the occurrence of short and recurring, 
but strongly super-critical (i.e. super-Eddington) accretion episodes 
at high-redshifts, with rates as large as 50-100 times the Eddington 
limit MEdd = LEddKec\ where: 


_ 4;rGM.c 
l-Edd =- 

Kj 


( 1 ) 


In the definition of the Eddington luminosity, G is the gravita¬ 
tional constant and kj is the Thomson opacity. The radiative ef¬ 
ficiency of the gas inflow depends on the accretion rate. Recall¬ 
ing the definition of the Eddington luminosity (Eq. 1), we define 
fEdd = Fl,IMEdd^ i-e- the accretion rate normalized to the Edding¬ 
ton value. It is expected that if matter is accreted at moderate rates 
0.01 "^fEdd ^ the inflowing material creates a radiatively effi¬ 
cient, geometrically thin and optically thick accretion disk (Shakura 
& Sunyaev 1973). In this case the radiative efficiency depends only 
on the black hole spin, and varies from ~ 6% for Schwarzschild 
black holes to ~ 32% for maximally rotating ones (Thome 1974). 
If, instead, accretion occurs super-critically (fEdd >1) the structure 
of the accretion disk is modified because of advection: the energy 
produced in the disk is carried inwards, in the black hole, rather 
than being radiated away (see, e.g., Abramowicz & Eragile 2013; 
Lasota 2015). The disk thickness increases and the disk becomes 
geometrically thick. The most common solution proposed for such 
accretion flows is the “slim disk" (Abramowicz et al. 1988; Paczyn- 
ski & Abramowicz 1982; Mineshige et al. 2000; Sadowski 2009, 
2011; McKinney et al. 2014), radiatively inefficient and with a thick 
geometric structure, in which photon trapping is significant inside 
the trapping radius Rpf. 


Rs M 
£ ^Edd 


( 2 ) 


where R, is the Schwarzschild radius. Only a fraction of the pho¬ 
tons produced by the viscous process inside the accretion flow is 
able to free stream out of the trapping radius, because the photon 
diffusion time exceeds the time scale for accretion. Consequently, 
the effective radiation pressure acting on the surrounding gas is 
decreased (see e.g. Begelman 1978; Ohsuga et al. 2002) and the 
luminosity is only mildly (e.g., logarithmically) dependent on the 
accretion rate. Alternatives of the slim disk solution exist, e.g the 
ZEro-BeRnoulli Accretion (ZEBRA Coughlin & Begelman 2014) 
and the ADiabatic Inflow-Outflow Solutions (ADIOS Blandford & 
Begelman 1999; Begelman 2012) models , which allow for a frac¬ 
tion of the inflowing mass to be lost during the accretion process. 

The present work investigates, both analytically and numer¬ 
ically, the growth of high-redshift (z ~ 10) black hole seeds. The 
growth process may be controlled by the amount of gas avail¬ 
able in the halo (feeding-dominated) or by the radiative feedback 
(feedback-limited). The growth is feeding-dominated if the radia¬ 
tive back-reaction of the black hole is negligible: the rapidity of the 
process is mainly determined by the gas accretion rate that the host 
halo can provide. The black hole growth is more efficient and rapid 
if the flow is feeding-dominated, assuming that a sufficient amount 
of gas is present in its host halo. 

We devise a very general analytic model which is able to pre¬ 
dict the growth efficiency from the physical properties of the system 


formed by the black hole seed and its host halo. Eurthermore, we 
employ a radiation-hydrodynamic code to follow the growth pro¬ 
cess from small (0.002 pc) to large scales (» Rb, where Rb is the 
Bondi radius), spanning a spatial dynamic range of four orders of 
magnitude, with special emphasis on the properties of the inner re¬ 
gions of the host halo, providing most of the accretion material to 
the central object. The growth is monitored as a function of differ¬ 
ent parameters, namely: the accretion model (radiatively efficient or 
inefficient), the density profile of the halo p(r) and the initial mass 
of the seed Mq = 10^“® Mq. 

This study is a natural follow-up of the paper by Pacucci 
& Eerrara (2015), where the authors have simulated in great de¬ 
tail the accretion process onto a z ~ 10 black hole seed of initial 
mass 10^ Mq, embedded in a dark matter halo with a gas mass 
of ~ 10^ Mq and with extreme density conditions (a number den¬ 
sity of hydrogen particles at ~ 0.1 pc from the center of the halo 
~ 10^ cm“^), finding that in ~ 142Myr about 90% of the gas mass 
of the halo has been accreted onto the compact object. 

In this general framework, we aim at clarifying several as¬ 
pects, including (i) do the radiatively inefficient accretion models 
provide an effective way to rapidly increase the mass of the seed? If 
so, why? (ii) What is the final fate of black hole seeds as a function 
of their initial mass? How long can they accrete? (iii) Is it possible 
to predict the growth efficiency of the accretion process from the 
physical properties of the (black hole -l- host halo) system? In the 
cosmological context, answering these questions will provide some 
insights into the mass that has been locked inside black holes during 
the cosmic evolution. This would be of great interest to understand 
both the formation of high-redshift SMBHs and their remnant pop¬ 
ulation. 

The outline of this paper is as follows. In §2 we briefly de¬ 
scribe the physics and equations of the radiation-hydrodynamic 
problem, along with the initial conditions for the density pro¬ 
files. In §3 we present our analytic model for the black hole 
growth, while in §4 we show the results of our simulations. Ei- 
nally, in §5 we provide some further discussion and the conclu¬ 
sions. Throughout, we adopt recent Planck cosmological param¬ 
eters (Planck Collaboration et al. 2015): (Q.„,,QA,^h^h,ns,crg) = 
(0.32,0.68,0.05,0.67,0.96,0.83). 


2 PHYSICAL AND NUMERICAL IMPLEMENTATION 

The present study employs a series of radiation-hydrodynamic sim¬ 
ulations to test the predictions of our growth model. Our code is de¬ 
signed to perform a fully consistent treatment of uni-dimensional 
spherically-symmetric hydrodynamic equations and a simplified, 
frequency-integrated, version of radiative transfer equations. While 
the spherical symmetry is an idealization of a real accretion flow, 
several works have shown that ID simulations provide a reliable 
description of many of its most important features. For instance, in 
Novak et al. (2011) the authors performed a comparative analysis 
between the outputs of a code run in ID and 2D, finding similar 
results in terms of black hole growth and duty cycle. The main 
difference, for which the multi-dimensional approach is a signifi¬ 
cant improvement, concerns the fact that the additional degrees of 
freedom allow classical instabilities (e.g. the Rayleigh-Taylor and 
the Kelvin-Helmholtz ones) to operate. Their net effect is to pro¬ 
duce a somewhat higher accretion rate, a less effective feedback 
and a more irregular pattern of bursts, compared with the ID case. 
Notwithstanding the dimensionality of the code, our implementa¬ 
tion cannot take into account the full complexity of the accretion 


MNRAS 000, 1-12 (0000) 




The Growth Efficiency of High-Redshift Black Holes 3 


flow, which only a more advanced treatment of angular momentum 
transport would allow. Due to the triaxiality of the host halo, the an¬ 
gular momentum field is extremely complex and variable in time, 
at every location inside the inflow (see Choi et al. 2013, 2015). 
Overall, the outward angular momentum transfer is very efficient 
(Choi et al. 2015), due to the gravitational torques induced by both 
the dark matter and the gas distributions. The gas loses its angular 
momentum efficiently and flows well beyond its centrifugal bar¬ 
rier: the gas that reaches the black hole is expected to have low 
angular momentum. Therefore, despite its simplifications, our ap¬ 
proach helps in acquiring physical insight on the process of black 
hole growth. In this Section we present a general overview of the 
most important aspects of the physical implementation, while the 
interested reader is referred to the paper Pacucci & Ferrara (2015) 
for a detailed description of the code. 

The domain of our simulations* spans approximately from 
0.002pc to 20pc. The characteristic spatial scale for accretion is 
the Bondi radius Rg: 


Rb 


GM. 

c^, , 

5 ( 00 ) 


(3) 


Here, Cj(oo) is the sound speed at large distances from the black 
hole, defined as: 


C.?(oo) — 



(4) 


where 7 = 5/3 is the ratio of specific heats, R is the gas constant. Too 
is the gas temperature at large distances and /i = 1.15 is the mean 
molecular weight for a primordial gas. For instance, the Bondi 
radius corresponding to a seed with initial mass M, = 10^ Mq 
with Cs(oo) ~ 12kms“* (i.e. Too ~ 10^K) is Rg = 3pc ~ 10“^f?vir ~ 
10^fJs> where is the virial radius of the halo. Our spatial 
range, spanning ~ 4 orders of magnitude, covers the entire range 
of Rg corresponding to an initial mass of the seed within the range 
10^-*= Mq. 

In the hydrodynamic module we solve the standard system 
of ideal, non-relativistic Euler’s equations (conservation laws for 
mass, momentum and energy, neglecting viscosity, thermal con¬ 
duction and magnetic fields) for a primordial (H-He) composition 
gas with helium fraction Yg = 0.24665 (Planck Collaboration et al. 
2015) and no metals, spherically accreting onto a central black 
hole, supposed at rest and already formed at the beginning of the 
runs, with a given initial mass Mq = M,(t = 0). The infalling gas 
has zero angular momentum with respect to the black hole. 

The forces acting on the gas are: (i) the thermal pressure, (ii) 
the gravitational pull of the black hole and (hi) the radiation pres¬ 
sure. The thermal pressure is given by the usual equation for ideal 
gases: 


pRT 


(5) 


where p is the mass density. The gravitational acceleration at the 


' We performed a series of convergence tests on the extension of the spatial 
range (see Pacucci & Ferrara 2015) which confirmed that the main outputs 
of our simulations (e.g. duty cycles and accretion rates) do not depend on it 
as long as (i) it covers a sufficiently large radial range around Rg (e.g. from 
~ O.lRfi to ~ 2Rb) and (ii) the centrifugal radius, where the accretion disk 
starts to form, is not resolved. 


distance r from the central object is: 

, , GM.fi) 
g(r,t) = -^—. 


( 6 ) 


The value of the black hole mass M.it) changes with time, due to 
gas accretion, with the following set of rules, where M, = 4;rr^p|v|: 


{M.{t)tf) « v(r = ro,0<0 

\M.(f) = Mo-Kl-6)/jM.dt 

where v is the velocity of the gas and e ranges^ from e = 0.057 
for a Schwarzschild (i.e., non-rotating) black hole to e = 0.32 for a 
maximally rotating object (see Thorne 1974). The mass flow M, is 
computed at the accretion radius ro (the innermost cell of our spa¬ 
tial grid) and it is equivalent, by definition, to the accretion rate onto 
the black hole. The acceleration caused by the radiation pressure is: 


^radd) ~ 


K(p,T)L{r) 

4i:Ec 


( 8 ) 


where L is the emitted luminosity and K{p, T) is the opacity of the 
gas, which includes the Thomson term (with the additional inclu¬ 
sion of a temperature dependence, as in Begelman et al. 2008) and 
bound-free terms. In our code, all radiation-related quantities are 
integrated over frequencies: matter and radiation are coupled via 
Thomson (electron) scattering and bound-free interactions. The ra¬ 
diative transfer employs a simplified two-streams approximation 
presented in Novak et al. (2012), to which the interested reader is 
deferred for a detailed description. 

The relation between the emitted luminosity Lq = Rdk) and 
the accretion rate M, is: 


Lo = ec^riM.) ifv(ro )<0 
To = 0 if v(ro) > 0 


where T'fM.) is a generic function of the accretion rate. In the sim¬ 
ple case of radiatively efficient accretion (treated in Pacucci & Fer¬ 
rara 2015) the following relation holds: 

7 -= fcdd- ( 10 ) 

CEdd 

In the radiatively inefficient accretion mode (for instance, the slim 
disk solution) the luminosity L depends on fg^hi as in the following 
prescriptions (Mineshige et al. 2000; Volonteri et al. 2014): 


L 

^Edd 


fsdd 

25 


ifEdd < 50); 


( 11 ) 



l-Edd 


1 -t-ln 


fEdd 

50 


(fEdd >50). 


( 12 ) 


In this scenario, only a fraction of the emitted luminosity escapes 
to infinity, and the effective radiation pressure depends only weakly 
on the accretion rate: note that L = ILgdd for fgdd = 50. In the ra¬ 
diatively inefficient accretion mode the matter-radiation conversion 
efficiencies are e = 0.04 for fg^id < 50 and: 


e 


1 / fEdd 'l * 

25^ 50 j 


1 -rln 


fEdd \ 

50 ] 


<0.04. 


(13) 


^ We assume that the matter-radiation conversion factor e is numerically 
equal to the accretion efficiency r], i.e. all the gas reaching the inner bound¬ 
ary is actually accreted by the black hole (see e.g. Sadowski et al. 2014; 
McKinney et al. 2014; Jiang et al. 2014). 
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for fEdd > 50. 

Due to radiation pressure, the accretion may occur in an inter¬ 
mittent manner. Therefore, it is useful to introduce the duty cycle, 
defined as the fraction of time spent accreting within a given time 
frame of duration t,ot'- 


D=l 


Udle 

hot 


(14) 


Here, is the idle time, i.e. the time spent without accretion tak¬ 
ing place. 


2.1 The initial density profile 


In this work we simulate the spherically-symmetric gas accretion 
onto a black hole seed at the center of a dark matter halo of virial 
temperature r,,,> ~ 10“* K and total mass (dark matter and baryons) 
Mh = 6.7 X 10^ Mq at z = 10. For r « as in our computational 
domain, most of the mass is baryonic. We assume that the gas ini¬ 
tially follows the isothermal (T ~ 10^ K) density profile derived 
from the simulations in Latif et al. (2013a), well approximated by 
the functional form: 


p{r) = 


PO 

1 -I- (rja)^ ’ 


(15) 


where a is the core radius. 

In order to study how the black hole growth depends on the 
host halo, we implemented two different density profiles, schemat¬ 
ically shown in Fig. 1. A high density profile (HDP) with a central 
density po = 10“*^ gcm“^ and a core radius a - 1.6x10“^pc and a 
low density profile (LDP) with a central density po = 10“** gcm“^ 
and a core radius a = 2pc. Both density profiles yield a total bary¬ 
onic mass Mgas ~ 10^ Mq over the entire spatial domain. In the 
HDP case the spatial domain has been slightly enlarged, with re¬ 
spect to the LDP one, to fulfill this condition. The yellow-shaded 
area in Fig. 1 shows that the LDP may be interpreted as the den¬ 
sity profile resulting after the formation of a black hole of mass 
~ 10^ Mq at the center of the halo (see Latif et al. 2013b, 2014). 
Consequently, the HDP may be interpreted as the density profile 
resulting after the formation of a very small (< 10* Mq) black hole 
seed, which leaves the matter distribution of the halo almost un¬ 
altered. In summary, we study four models: HDP with radiatively 
efficient accretion, HDP with radiatively inefficient accretion, LDP 
with radiatively efficient accretion, and LDP with radiatively inef¬ 
ficient accretion. 


3 ANALYTICAL INSIGHTS 

In this Section we present our analytic model to describe the growth 
of high-redshift black holes. We predict whether the mass inflow is 
feeding-dominated or feedback-limited, as a function of the main 
physical properties of the (black hole -l- host halo) system. This 
analytic framework is then employed to determine a time-evolving 
spatial scale outside of which the radiative feedback is highly effec¬ 
tive and outflows dominate. Moreover, we envisage the existence 
of a mass scale above which the accretion flow is always feeding- 
dominated. 


3.1 Modeling the mass inflow 

Our analytic model takes into account the following properties: 

(i) the density of the gas in the inner sections of the halo po, 

(ii) the mass of the black hole seed M,, (iii) the matter-radiation 



Figure 1. The two density profiles employed, both yielding a total baryonic 
mass of ~ 10^ Mq over the integration range. The LDP has a much larger 
core radius, but the central density is smaller by ~ 6 orders of magnitude. 
The yellow-shaded region contains a total mass of 10^ Mq and may be inter¬ 
preted as the gas which has been extracted from the HDP to produce a hlack 
hole of the same mass. Consequently, the HDP profile may be interpreted 
as the density profile resulting after the formation of a small hlack hole of 
mass M, < IO^Mq. 


conversion factor e, (iv) the Eddington factor fEdd- Starting from 
quasi-static conditions, the accretion rate increases, due to the black 
hole’s gravitational pull. The build-up of the accretion rate causes 
an increasing emission of radiation, which eventually might be able 
to stop the gas inflow, or even invert its velocity: this depends on 
the intensity of the radiation field and on the inertia of the gas. 
The rate at which the linear momentum of the gas is changed relies 
on the comparison between two time scales: (i) the feedback time 
scale, tfh, defined as the time needed by the radiation pressure to 
significantly (i.e. by a factor e) change M,, and (ii) the accretion 
time scale, tacc, which estimates the time needed to consume the 
gas mass inside the inner regions of the accretion flow. 

Given the general expression for the acceleration caused by 
the radiation pressure (Eq. 8) and the relations between the lumi¬ 
nosity L and the Eddington factor fEdd for radiatively efficient and 
inefficient flows (Eqs. 10, 11, 12), we calculate tfh. The impulse 
theorem (radial component) states that dp = Fdt, where p is the 
linear momentum per unit volume of the gas, F is the force applied 
on the gas per unit volume and t is the time: 

K L 

d(pv) = pUraddt =p- - fdt. (16) 

Anc 

This equation needs to be evaluated at some radius representative 
of the accretion process onto the black hole. In a simulation, this 
radius is the innermost cell of the spatial grid, where the gas is 
assumed to be accreted by the black hole. From a purely theoretical 
point of view, it may be considered as the radius where the accretion 
disk forms. In either case, we designate this spatial scale rg, the 
accretion radius. Evaluating the previous equation at r = rg gives: 

ii(4;rrQpovo) = Lodt. (17) 

^ c 

Here, Lq is computed with Eq. 9 and the mass flux is computed 
at the accretion boundary, corresponding to the accretion rate M,. 
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Therefore: 

dM.^—Lodt. (18) 

c 

We solve this ordinary differential equation with the initial con¬ 
dition M,(t = 0), assuming that the accretion rate decreases with 
time: in the case of an increasing mass flux, only the sign of the 
exponential solution changes. The general solution is: 

M.(t) = M.(t = Q)e^’l‘f>’, (19) 


where 


and 


d 

epoKQC ’ 


lA = fsdd 


t-Edd 

L 


( 20 ) 


( 21 ) 


In the standard radiatively efficient scenario e = 0.1 and lA = 1, while 
in the radiatively inefficient case e = 0.04 and if/= 25 (for fEdd < 50) 
and e < 0.04 and if/>25 (for fEdd ^ 50). The feedback time scale is 
not related to the black hole mass M,, but only to the mass of the 
inflowing gas, via po. In the slim disk case, for a given density po, 
the feedback time scale is > 60 times longer: the system reacts to a 
modification of the accretion rate in a much slower way. 

To calculate tacc, we estimate the time scale for the consump¬ 
tion (due to accretion) of the gas mass inside some radius r, given 
an accretion rate M ,: 


, 22 ) 

M. iira/EM M. • 

where on the right-hand side we have parametrized the luminos¬ 
ity through the Eddington factor fEdd- The accretion time scale^ 
is inversely proportional to the black hole mass M, and directly 
proportional to the gas mass within the radius r. 

Defining the ratio of the two time scales as 'r(r,t) = tf^ftacc^ 
a transition radius, rj, exists such that taccdr) = tftivT)- The 
flow is feeding-dominated in the region where r <K rj, tacc ^ 
T{r,t) » 1 and gas is easily available for the accretion. The flow is 
instead feedback-limited where r»rj, tacc 3> f/;,, 7”(r,r) « 1 be¬ 
cause outflows dominate over inflows and accretion proceeds in an 
intermittent way, or can be even halted if the gas reservoir is empty. 
The parameter 7”(r, t) allows us to determine whether an accretion 
flow is feeding-dominated or feedback-limited. The expressions for 
rj and T(r,t) are as follows: 


and 


n 


. 'iGfEddM, 

r - 

(epQKQCY 


T{r,t) = if/ 


^GfEddM. 

(epQKQcYd 



(23) 


(24) 


Note that the black hole mass M, increases with time, while po 
decreases as gas is consumed, therefore the transition radius in¬ 
creases with time: a larger fraction of the host halo enters the 
feeding-dominated region, and accretion becomes progressively 
easier. Moreover, we define T{rQ,t) = Tq. 

Eq. 24 shows that the efficiency of an accretion flow depends 
on several variables that we discuss in turn. Eirstly, the smaller 


^ The accretion time scale is also very well approximated by a fraction 1/3 
of the crossing time for a particle at some radius r. MfM = r/{3r). 



Figure 2. Comparison, at t = 0, between the accretion time scale (tacc) and 
the feedback time scale (tfi,) for the accretion flow onto a black hole of 
initial mass 10^ Mo in the radiatively efficient and LDP case. The black line 
is the value of T” (r,to)- The yellow-shaded area is the region of integration 
in the present work, while the orange-shaded region indicates the transition 
region, where 0.1 < 7“ <10. The inner and outer boundaries, the Bondi 
radius and the virial radius of the halo are shown. The black star indicates 
the position of the transition radius. 


the internal density of the halo, po, the longer is tfh- This is be¬ 
cause when the gas density is low the physical accretion rate on 
the black hole is small, and, at a given black hole mass, radia¬ 
tion pressure is less effective. It will be easier to have a feeding- 
dominated flow in the LDP than in the HDP case. Secondly, a 
smaller radiative efficiency e yields a less effective radiation pres¬ 
sure, and consequently radiatively inefficient accretion would be 
more likely feeding-dominated. Thirdly, a larger Eddington factor 
fcdd oc V implies a higher inward linear momentum of the gas. 
In this case, the inward velocity of the gas mass is less easily in¬ 
verted. Lastly, the black hole mass M,: since MEdd a smaller 
mass corresponds to a smaller physical critical accretion rate: in a 
given halo a small black hole is more likely to be fed at rates that, 
for its mass, give rise to high radiation pressure. A small black hole 
is therefore more likely to find itself in the feedback-limited regime. 

As an example, Eig. 2 shows a comparison between the ac¬ 
cretion time scale (tacc) and the feedback time scale (tfi,), both 
computed at f = 0, for the accretion flow onto a black hole of ini¬ 
tial mass 10^ Mo embedded in an LDP in the radiatively efficient 
case (Eq. 10). The transition radius rj ~ 2x 10“^pc is larger than 
the accretion radius (ro ~ 2 x 10“^ pc), so we expect the black hole 
growth to be feeding-dominated. Moreover, an overall modification 
of the spatial velocity profile of the accretion flow should be visible 
around the position of the transition radius, in the spatial range that 
we schematically call transition region, where 0.1 < 7” < 10. These 
effects are discussed in Section 4, dedicated to the results of our 
numerical simulations. 

In the HDP case the transition radius would be smaller than the 
accretion radius and the growth would be feedback-limited. On the 
contrary, in the radiatively inefficient cases, both HDP and LDP, the 
transition radius would be even larger than the one shown in Pig. 2, 
leading to a more extended feeding-dominated region. 
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3.2 Relation between 2) and T 3 

In the following we derive an analytic relation between our model, 
through the quantity Tq (Eq. 24), and D, the duty cycle for the 
black hole growth (Eq. 14), which is a phenomenological way, 
computable only a posteriori, to describe if the gas inflow is contin¬ 
uous. Even Adaptive Mesh Refinement cosmological simulations 
cannot resolve the typical spatial scale of accretion, therefore they 
have to resort to some kind of sub-grid prescriptions to model the 
black hole growth, like assuming that they continuously accrete at 
the Bondi-Hoyle-Lyttleton rate (Bondi 1952), capped at the Ed¬ 
dington rate (see e.g. Springel et al. 2005; Di Matteo et al. 2008; 
Dubois et al. 2013; Costa et al. 2014; Dubois et al. 2015). Our 
model may provide in such cases more realistic values for the duty 
cycles of these sources. 

Calling tijig the fraction of t,o, during which the black hole is 
not accreting, it is possible to show that this relation holds: 

— = g-'/f^ (25) 

hot 

by solving the following differential equation (where ^ is the time 
scale for the variation of M,): 

-^ = ^. (26) 

dt 'l.tacc 

The factor 2 accounts for the fact that, once the radiation pressure 
exerts an action on the infalling gas for some time At, a time 2Af is 
needed to re-establish the accretion rate preceding the acceleration 
(see Eq. 19). We obtain: 

D = (27) 

With our definition, DiT 3 = 1) ~ 0.4. A proof of the validity of Eq. 
27 is provided in Fig. 7, described in section 4.2.1. 


3.3 The black hole - host halo connection 


The transition radius, which separates the feeding-dominated re¬ 
gion (r « rj) from the feedback-limited region (r » rj), increases 
with time (cf. Eq. 23). Consequently, if the black hole growth is 
feeding-dominated at rg and f = 0, it will always be so. Asking that 
W(ro, to) > 1 translates into a black hole mass above which the flow 
will always be feeding-dominated: 


Mrrit = - 


3x 10^ 
^fsdd 
ro 


1 f V 

7 Po 'i" 

\0.1/ 

\5x 10-^®gcm-3/ 


X 


10-4 pc 


Mp 


(28) 


Here, we considered ro ~ 10 “^ _ pc ~ 20 - 200AU as a typical 
spatial scale for accretion disks. 

In the HDP case for the standard accretion scenario the previ¬ 
ous limit reads: 


T(ro,to) > 1 if M. > Merit ~ 3 X 10® Mo (HDP - Std), (29) 


while both in the LDP case and in the slim disk accretion scenario 
(both density profiles) the limit is negligible: 

T~FoJo)> 1 if iff* > iffcriY ~ 10- IOOMq (othercases). (30) 


The physical meaning of Merit requires a clarification, since 
one normally expects that feedback halts the black hole growth 
above a given mass, rather than below (see e.g. Silk & Rees 1998; 
King 2003, 2010). The meaning of the lower limit we find is that 
when M, > Merit the accretion rate needed to exert a sufficiently 


strong feedback is so high, for the halo in question, that the ac¬ 
cretion flow cannot produce it. In other words, the inflow rate is 
determined by the halo properties, and is, at least initially, indepen¬ 
dent of the black hole mass. If a given halo provides the same M, to 
a small black hole or a large black hole, it will be the smaller black 
hole that will reach the Eddington limit first, having its growth 
stunted. As a general result, smaller black hole seeds should en¬ 
counter great challenges during the first stages of the growth, char¬ 
acterized by outflows and very low values of the duty cycle. This 
effect could play an important role at high redshifts, where black 
hole seeds of different mass may form from the same host halo, de¬ 
pending on the thermal and radiative properties of the environment. 


4 NUMERICAL RESULTS 

The following two subsections describe our numerical simulations 
and their analysis through the growth model outlined so far. 


4.1 Inflows and outflows 

As an example, we discuss a simulation of the accretion flow onto 
a black hole of initial mass 10 ® Mq in the radiatively efficient and 
LDP case. Starting from static and isothermal conditions, the grav¬ 
ity of the black hole rapidly pulls in the gas, building up the ac¬ 
cretion flow. Fig. 3 shows the time evolution of the velocity and 
temperature spatial profiles, which may be interpreted as the dy¬ 
namical and thermal counterparts of Fig. 2, bearing in mind that 
rrit) always increases. The transition and the outflow regions are 
defined as the set of radii r where, at least once during the time evo¬ 
lution of the system, the relations 0.1 <T(r) < 10 and T(r) <0.1 
respectively hold. In the region close to the accretion boundary the 
inflow is very smooth, with inward velocities up to ~ 15kms-* 
and temperatures rising up to ~ 2.4 x lO^K. In the transition re¬ 
gion, the flow starts to be disturbed by radiative feedback, which 
becomes more effective due to the increase of tacc, with frequent 
velocity inversions and a more complex temperature profile. In the 
outflow region the radiative feedback is dominant, with large (up to 
~ 5kms-*) outflowing velocities and a temperature profile which 
reconnects to the thermal floor (T ~ 10^ K) of the host halo. 

For the same simulation, Fig. 4 shows the time evolution of the 
mass flux at the inner boundary (i.e. the accretion rate), at the outer 
boundary and inside the transition region, in which the computed 
values are a spatial average over the cells belonging to this region 
at each time of the simulation. The inner region is characterized by 
a slow, but constant, accretion of order ~ 2.0x lO-^Moyr'*. The 
outer region is swept by large outflows, whose magnitude increases 
with time, reaching a peak of ~ 0 . 8 Moyr-*, while the transition 
region is characterized by a mild outflow. 

In the HDP case (for the same accretion scenario), rj <ro and 
the growth is feedback-limited: the accretion is discontinuous and 
the transition region would extend down to the accretion boundary, 
while the outflows would be more intense. On the contrary, in the 
radiatively inefficient cases, both HDP and LDP, rj » ro: the accre¬ 
tion flow would be smooth and continuous over a large fraction of 
the spatial domain, with the transition region starting to be visible 
only close to the outer radius. 


4.2 The growth efficiency 

In the following we numerically test the predictions of our model 
with 16 runs, in order to adequately explore the range of density 
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Figure 3. Velocity and temperature spatial profiles for the accretion flow 
onto a Mq = 10^ Mo black hole embedded in a LDP density profile, accret¬ 
ing in the radiatively efficient mode. The total integration time is = 10^ yr 
(ti„i = iAt with At = 10“* yr and i = 0, 2,4, 6, 8, 10). The initial conditions are 
marked with the label “IC". The transition and outflow regions are defined 
by 0.1 '£'?”('■)< 10 and 7”(r) < 0.1 respectively. The inner and outer bound¬ 
aries are shown with orange dashed lines. 

profiles (HDP and LDP), black hole mass seeds (IO^^^Mq) and 
radiative efficiencies (L oc and L oc In/gj^). The simulations 
are run for ~ 10^ yr to allow all the simulated accretion flows to 
reach steady-state conditions without arriving at a complete deple¬ 
tion of the gas reservoir. The parameters reported in the figures, 
{T~{r(i)), ifEdd) and {O), are an average over the entire simulation 
time, while the parameter AMIMq expresses the total mass growth 
at the end of the simulation, normalized to the initial value of the 
seed mass. Fig. 5 shows the results in the radiatively efficient case, 
while Fig. 6 refers to the radiatively inefficient one. The top panels 
refer to HDP, while the bottom ones to LDP. 

4.2.1 The M. = 10^^® Mq runs 

The HDP case in the radiatively efficient scenario (Fig. 5, top panel) 
is particularly interesting because 7~ crosses unity twice, allowing 
us to test the O — 'T relation. The simulated masses are all below 



Time [10“ yr] 

Figure 4. Time evolution of the mass flux computed at the inner boundary, 
outer boundary and inside the transition region highlighted in Fig. 3. There 
is a slow but constant (~ 2.0x 10“^ Moyr^') accretion at the inner boundary, 
while the outer region is characterized by a strong outflow, with a mass flux 
up to ~ 0.8 Mq yr“*. The transition region is characterized by a mild outflow, 
on average. 


Merit ~ 3 X 10® characteristic of this density profile and radiative ef¬ 
ficiency for fEdd = 1 . therefore one expects the flow not be feeding- 
dominated at all times. However, Pacucci & Ferrara 2015 show that 
super-critical accretion is feasible, for short times, also in radia¬ 
tively efficient scenarios. Indeed, in the mass range IO^^'^Mq the 
entire galaxy is in outflow {rj ~ fo), but the radiation pressure is 
active only when the gas accretion is underway, while during the 
remaining time the inflow huilds up again. As a consequence, the 
radiation pressure is unable to sweep all the gas away from the ac¬ 
cretion boundary, and a small physical accretion rate is sufficient 
to grow such small black holes, so the flow can sustain, on aver- 
age, fEdd ~ 25, with duty cycles 2) ~ 0.8 - 0.9. The relevant Merit, 
therefore, is not the one for f£dd ~ 1. but for f^dd ~ 25. The val¬ 
ues of T and rj slowly increase as the initial mass of the seed M, 
approaches Merit ~ L2x 10® Mq for f^dd ~ 25 (see the left-most 
vertical line in Fig. 5, top panel). 

In the range 10®~®Mo a super-Eddington flow is no more sus¬ 
tainable, since the radiation pressure is progressively more power¬ 
ful: the flow stabilizes at f^dd ~ 1- For an Eddington-limited flow. 
Merit ~ 3 X 10® Mo (see the right-most vertical line in Fig. 5, top 
panel). In the process of going from fEdd ~ 25 to fEdd ~ 1 the 
flow becomes mildly feedback-limited {T ~ 0.7 for M, = 10® Mo, 
2) ~ 0.3 -0.5). The physical accretion rates in the HDP case are 
within the range 10“® - lO^'^Moyr^* (see Fig. 5 for further de¬ 
tails). 

In the LDP case of the radiatively efficient scenario, the crit¬ 
ical mass value Merit is < IO^Mq: as a consequence the growth 
is feeding-dominated ( 2 ) ~ 1 and 7” > 10 ^) for all runs and accre¬ 
tion rates are stable, close to the Eddington value {fEdd ~ !)• The 
physical accretion rates are within the range 10 “^- 10 “®Moyr“', 
increasing with the mass M, just as the Eddington rate. 

In the HDP case of the radiatively inefficient scenario (Fig. 
6 , top panel) the situation is similar, because Merit < 10 ® Mq: so 
tt » rg and the growth is always feeding-dominated (fT up to ~ 
10^). Accretion is not capped at the Eddington rate, so we reach 
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Seed Mass [M^] 



Seed Mass [M^] 



Seed Mass [M^] 


Figure 5. Radiatively efficient case - Comparison between the parameters 
used to describe the accretion flow (7~, AM/Mq, fEdd and £)), for black 
hole seed masses in the range 10^“^ M© and two density profiles (HDP 
above and LDP below). The vertical purple lines show the values of Merit 
for fEdd - 25 and fEdd = from left to right. The physical accretion rates 
are, in ascending order of mass: 4.Ox lO^'^Moyr^^, 4.5 x lO“^M 0 yr“^, 
1.5 X 10^^ Mo yr^', 6.5 X 10^^ Mq yF for the HDP, and 2.0 X 10^’ Mq yr“‘ , 
2.0x lO^'^Moyr-*, 2.0x IQ-^Moyr-*, 2.0x IQ-^Moyr^* for the LDP. The 
simulation with Mq = 10^ Mo shows a very large value for the physical ac¬ 
cretion rate, ~ 2.0 X 10“* Mo yr“*. See the main text for further details. 


rates as high as fEdd ~ 300, leading to a super-Eddington emitted 
luminosity L ~ ShEdd (Eq- 12). The physical accretion rates reach 
large values, with a peak of 14Mo yr“*. 

Similarly, in the LDP case of the radiatively inefficient sce¬ 
nario (Fig. 6, bottom panel) all accretion flows are feeding- 
dominated, with physical accretion rates within the range 10“*- 
10“^Moyr“*. While the Eddington factor is fEdd ~ 5 for all 
masses, the emitted luminosity is strongly sub-Eddington, due to 
the low mass density of the LDP scenario, which is unable to pro¬ 
duce a mass inflow leading to super-critical luminosities. 

In summary, radiatively inefficient accretion allows for largely 
feeding-dominated growths of the central black hole, while stan¬ 
dard accretion scenarios may be feedback-limited when the gas 
density is very high. Moreover, the main difference between the 
two density profiles, notwithstanding the accretion mode, is related 


Figure 6. Radiatively inefficient case (slim disk) - Same as Fig. 5. The 
physical accretion rates are, in ascending order of mass: 7.0 X 10“^ Mq yr“* , 
7.0x lO-^Moyi-*, 1.5Moyr-‘, I.4x lO'Moyr-* for the HDP, and 2.0x 
10“''Moyr-‘, 2.0 X Moyr^*, 2.0 X lO^^Moyr^*, 1.5 X lO^* Moyr^' 
for the LDP. See the main text for further details. 


to the accretion rates that they can sustain to feed the black hole: 
the HDP produces an accretion rate up to fEdd ~ 300, decreasing 
dramatically the growth time of the central object. 

With the average values of T and D computed so far for the 
16 runs, we are in the conditions of testing our model for the D — 
T relation, given by Eq. 27. Eig. 7 shows that our analytic model 
offers a very good fit to the points in the (2), 7”) plane. 

In order to have a general overview of the mass accreted 
(across the inner boundary) and ejected (across the outer boundary) 
in the different simulations performed. Fig. 8 provides a mass bal¬ 
ance summary. The mass growth (AM, h M,(t) — Mq) for a black 
hole seed of initial mass 10^ Mq is shown with solid lines, while the 
dashed lines indicate the mass ejected through outflows. The black 
solid diamonds indicate the final mass that the same seed would 
reach accreting continuously at the Eddington rate (2) = 1, fEdd =1) 
with e = 0.1. The ejected mass is larger than the accreted mass, ex¬ 
cept for the HDP in the slim disk scenario (brown line), due to the 
large values of fEdd available under these physical conditions. The 
mass growth in the standard radiatively efficient scenario is smaller 
than or equal to the one predicted by a continuous accretion at the 
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Figure 7. Data points in the (S,7”) plane, obtained from the 16 runs per¬ 
formed (Figs. 5 and 6). The overplotted line is the theoretical relation be¬ 
tween D and 7~, discussed in the main text (Eq. 27). 


Eddington rate, while in the slim disk scenario it is 2 - 3 orders of 
magnitude larger (in particular, the slim disk and HDP accretion 
provides a final mass which is ~ 1000 times larger). 

In a preceding work (Pacucci & Ferrara 2015), the authors 
found that only ~ 10% of the halo mass is expelled with outflows 
during the accretion process. In the present work, during the first 
10^ yr of evolution of a 10^ Mq seed embedded in a HDP halo the 
outflowing gas is ^ 85% of the total. This difference is explained by 
three reasons. Firstly, the relative importance of outflows decreases 
with time, due to the fact that the gradient of rj-fr) is positive; the 
simulation presented in Pacucci & Ferrara (2015) is ~ 1000 times 
more extended in time than the ones analyzed here. Secondly, the 
spatial range of the present work is much larger, then more apt to 
probe the outflow regions. Lastly, the halo density profiles are dif¬ 
ferent, since the one employed in Pacucci & Ferrara (2015) had a 
higher central density. Most of the baryonic mass of a r,,,y ~ 10"^K 
halo at z ~ 10 was confined inside a sphere of radius ~ 3 pc, lead¬ 
ing to values of the optical depths of order Nh ~ 9 x 10^^ cm~^ at 
the beginning of the collapse, which may have strongly reduced the 
impact of outflows. 


4.2.2 The M, = 10^ Mq run: a test case for Merit 

In the radiatively efficient case with HDP (Fig. 5, top panel) the 
value of 7” reaches a minimum (~ 0.7) around M, ~ 10^ Mq and 
rises up to ~ 60 again for M, ~ 10^ Mq, allowing us to test the 
value of Merit ~ 3 x 10^ Mq (Eq. 28) predicted for this scenario. 
For M, Merit, the Eddington rate is so large (> 0.3Moyr“*) that 
the accretion flow from the halo simply cannot provide it (the free- 
fall rate at f = 0 is ^ 0.2Moyr“*). As a consequence, the Eddington 
factor remains below unity (fEdd = 0.7) while O ~ 1. In this case, 
the physical accretion rates are very large (~ 0.2Moyr“*), a factor 
~ 100 higher than the values for other seed masses (see the cap¬ 
tion of Fig. 5). This feeding-dominated simulation has produced 
continuous accretion with very large rates and with sub-Eddington 
luminosities. 




Figure 8. Mass accreted (solid lines) and ejected (dashed lines) for a black 
hole seed with initial mass 10^ Me, embedded in a HDP (top panel) and in 
a LDP (bottom panel), in radiatively efficient and inefficient scenarios. The 
black solid diamonds indicate the final mass that the same seed would reach 
accreting continuously at the Eddington rate with e = 0.1. 


4.3 The final black hole mass 

The integration time for all simulations is much shorter than the 
typical evolutionary time scale for these systems (Pacucci & Fer¬ 
rara 2015; Pacucci et al. 2015). In this section we estimate the de¬ 
pletion time tend needed to void the inner regions of the host halo 
from its gas content, due to both gas accreted by the black hole, and 
ejected through outflows. Furthermore, from tend we can provide a 
rough estimate of the final black hole mass, by extrapolating the 
average accretion rates up to tend- 

Extrapolating the lines in Fig. 8, we find the depletion time 
tend with the following condition: 

March — hnd) + Mgjit = tend) — Mgas , (31) 

where Maech) is the accreted mass, Mgjit) is the ejected mass and 
Mgas = 10 ^ Mq is the total gas mass within our computational do¬ 
main. Table 1 provides a general outline of the accretion history 
for a black hole seed with initial mass M, = 10^ Mq in the four 
scenarios investigated so far, including the depletion time tend, the 
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extrapolated final mass of the black hole and its ratio with 

the initial baryonic mass of the host halo M,(tf,„d)IMgas- 

While in standard accretion scenario the typical time scale is 
~ lOOMyr (Pacucci & Ferrara 2015) and the black hole can ac¬ 
crete ~ 5% - 15% of the baryonic mass of the host halo, in radia- 
tively inefficient modes the growth is much more rapid and effi¬ 
cient. Specifically, in the slim disk scenario we predict that the evo¬ 
lutionary time scale is of order ~ lOMyr (Pacucci et al. 2015), with 
outflows playing a negligible role: the black hole is able to accrete 
~ 80% - 100 % of the host halo gas. 


4.4 A bimodal evolution of the black hole seeds 


Several works in literature (e.g. Silk & Rees 1998; King 2003, 
2010) have investigated various forms of the so-called M, -cr rela¬ 
tion, which provides an upper limit for the black hole mass: a com¬ 
pact object embedded in a halo with velocity dispersion cr ~ Vesc 
(the halo escape speed) can grow up to a mass given by the M, - cr 
relation, while the remaining gas is dispersed by radiation-driven 
(Silk & Rees 1998) and/or momentum-driven (King 2003) out¬ 
flows. The usual assumption adopted in these works is that the mo¬ 
mentum flux (MoufV) of the outflowing gas is comparable to the 
one in the Eddington-limited radiation field: Moutv ~ LEddF- In 
this work we have taken an alternative view where, comparing the 
time scales for gas infall and gas ejection as a function of radius, 
we prove that the momentum flux may be very different from the 
Eddington value: for instance, in the slim disk model where super¬ 
critical accretion rates may be associated with sub-Eddington lumi¬ 
nosities. 

Our approach predicts the existence of a critical black hole 
mass Merit above which the accretion is negligibly affected by out¬ 
flows: this, in turn, may lead to a bimodal evolution of the initial 
mass function of high-redshift black hole seeds. The lower-mass 
seeds (M, < Merit) would go through a feedback-limited growth, 
with recurring episodes of strong outflows which deplete the in¬ 
ner regions of the host halo from its mass content: the black hole 
cannot accrete more than a few percent of the gas reservoir. On 
the contrary, higher-mass seeds (M, > Merit) would go through a 
feeding-dominated growth, with outflows playing a negligible role: 
the black hole grows in mass very rapidly, possibly even consum¬ 
ing most of the host halo mass, reaching the SMBH stage early in 
time. For the LDP case, and for slim disk accretion in either density 
profile Merit is very low (~ 10- 100 Mq), and therefore of relevance 
only if black hole seeds are stellar-mass or so. For the HDP case, 
and standard accretion. Merit '-3x10^ Mq is relevant for most seed 
masses proposed in the literature. 

As a proof-of-concept of this bimodal development. Fig. 9 
shows the cosmological evolution, between z = 10 and z = 7 (the 
epoch when the first SMBHs are observed), of two initial mass 
functions for high-redshift black hole seeds: a simple flat distri¬ 
bution in the mass range Logio(Af,[M 0 ]) = 4.5-5.5 (top panel) 
and the initial mass function modeled (at z = 10) in Ferrara et al. 
(2014) (bottom panel). This evolution, far from being a precise pre¬ 
diction of the actual black hole growth, is a proof-of-concept based 
on the theoretical framework described in the present work. The 
basic equation for the mass growth is the following one: 


M.(f) = M.(? = 0 )exp 


DfEdd 


t 

0.045 Gyr 


(32) 


where the values for I)(M,) and fEdd(M,) are interpolated for each 
M, from the solid lines in Fig. 5, in the HDP case. The bimodal evo¬ 
lution is evident and the mass gap at Log\Q{M,) ~ 5.7 is expected 


0.07 


Flat initial distribution of high-z black holes 



0.20 


Realistic IMF for high-z black holes 



LogwiM,) 


Figure 9. Proof-of-concept bimodal evolution, between z = 10 and z = 7, of 
two initial mass functions for black hole seeds: a fiat one (top) and a more 
realistic one (bottom). The evolution is computed along the theoretical lines 
described in this paper, with values for T)(M,) and fedd^M.) interpolated 
from the HDP simulations in the standard accretion scenario. 


to rapidly spread during the cosmic time. Importantly, this effect 
does not depend on the shape of the initial mass function. 

In the HDP case with a radiatively inefficient accretion the 
bimodal evolution is expected to occur as well since, while O ~ 1 
for all masses, the value of fEdd does show a mass dependence (see 
the upper panel of Fig 6 ). 

In halos with a LDP density profile, in any accretion scenario, 
the bimodal evolution is not expected to occur, since fEdd and O 
are nearly independent of the seed mass (see bottom panels of Figs. 
5 and 6 ). Nonetheless, we expect this evolutionary effect to play a 
remarkable role in the growth process, since we believe that dark 
matter halos with a HDP density profile harbored at their center the 
black hole seeds with smaller masses {M, < 10^“^ Mq), the ones 
which are affected the most by feedback-limited growth. 
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Table 1. Analytic estimate of the accretion history for a black hole seed with initial mass M, = 10^ Mq in the four scenarios investigated so far. 


Accretion scenario [Myr] M.(4„rf) [Mq] 


Standard accretion - HDP 

225 

7.4x10^ 

1% 

Standard accretion - LDP 

110 

1.5x10'^ 

15% 

Slim Disk accretion - HDP 

14 

1 .0x10'' 

100 % 

Slim Disk accretion - LDP 

9 

8.3x10'’ 

83% 


5 DISCUSSION AND CONCLUSIONS 

The aim of this work is to provide a theoretical framework, sup¬ 
ported by numerical simulations, to describe the growth of high- 
redshift (z ~ 10) black hole seeds. The growth can be either feeding- 
dominated or feedback-limited. It is feeding-dominated if the radia¬ 
tive back-reaction of the black hole is negligible: the rapidity of the 
process is mainly determined by the gas accretion rate that the host 
halo can provide. In a feeding-dominated accretion flow the black 
hole: (i) increases its mass during most of its time evolution, and (ii) 
exerts a relatively small mechanical and/or radiative feedback on its 
surrounding gas reservoir. For these reasons, a feeding-dominated 
accretion is the desirable way to grow the black hole efficiently. 

Within this theoretical framework, we investigated, with the 
aid of ID radiation-hydrodynamic simulations, the growth of black 
hole seeds with a set of astrophysically-motivated initial conditions 
to explore the consequences of our model: seed masses in the range 
10^“® Mq, embedded in a dark matter halo of total mass (dark mat¬ 
ter and baryons) Mf, = 6.7 x 10^ Mq with two different density pro¬ 
files and with different prescriptions for the accretion efficiency, 
namely a radiatively efficient mode (e = 0.1) and a slim disk mode 
(e<0.04). 

Three points in particular are worthy of being retained from 
this work. 

• We confirmed that radiatively inefficient accretion modes (for 
instance the slim disk model) may ensure a continuous growth with 
rates largely exceeding the Eddington limit (reaching ~ 300MEdd 
in our simulations). Radiatively inefficient accretion flows allow for 
feeding-dominated growths of the central black hole, while stan¬ 
dard accretion scenarios may be feedback-limited with high values 
of the host halo gas density. The feedback time scale for these radia¬ 
tively inefficient modes is >60 times longer than in the standard ac¬ 
cretion scenario: the system reacts to a modification of the accretion 
rate in a much slower way because at a given accretion rate the pro¬ 
duction of radiation is reduced, decreasing the feedback effective¬ 
ness. In addition, we numerically proved the feasibility of accretion 
flows with sub-Eddington luminosities and super-Eddington rates. 

• We theoretically derived the existence of a time-evolving tran¬ 
sition radius, rj, which discriminates between feeding-dominated 
and feedback-limited growths. The transition radius, in addition, 
determines the spatial scale at which outflows take place, and pro¬ 
vides a mass scale. Merit, above which the black hole growth is 
always feeding-dominated. The critical black hole mass is 10 - 
10® Mq, depending on the accretion scenario and on the host halo 
properties. Consequently, we foresee the possibility of a bimodal 
evolution of the population of black hole seeds: low-mass seeds 
grow much less efficiently than high-mass ones. 

• Our model may be employed in modeling the growth of high- 
redshift black holes in large cosmological simulations, which can¬ 
not resolve the typical spatial scales of accretion. 

To conclude, our model’s aim is to study high-redshift accre¬ 
tion flows leading to the growth of the first black holes. In par¬ 


ticular, we provided the theoretical framework needed to under¬ 
stand why radiatively inefficient accretion models are likely to be 
a crucial ingredient in explaining the presence of SMBHs of mass 
~ 10®“*®Mo less than IGyr after the Big Bang (Mortlock et al. 
2011; Wu et al. 2015). In order to build such extremely massive 
objects already at z ~ 6.5, an Eddington-capped accretion would 
require, starting from a stellar mass {M, ~ 50Mo) seed, a constant 
Eddington flow. However, a low-mass seed would be hindered in 
its initial growth by its own feedback, making continuous accretion 
at Eddington levels unlikely (see also Johnson & Bromm 2007; Al¬ 
varez et al. 2009; Milosavljevic et al. 2009; Park & Ricotti 2012). 


6 ACKNOWLEDGEMENTS 

We thank Bin Yue for useful comments and suggestions. MV ac¬ 
knowledges support from a Marie Curie FP7-Reintegration-Grant 
(PCIGlO-GA-2011-303609). 


REEERENCES 

Abramowicz M. A., Fragile P. C., 2013, Living Reviews in Relativity, 16, 1 
Abramowicz M. A., Czerny B., Lasota J. R, Szuszkiewicz E., 1988, ApJ, 
332, 646 

Alexander T., Natarajan R, 2014, Science, 345, 1330 

Alvarez M. A., Wise J. H., Abel T., 2009, ApJ, 701, L133 

Begelman M. C., 1978, MNRAS, 184, 53 

Begelman M. C., 2012, MNRAS, 420, 2912 

Begelman M. C., Volonteri M., Rees M. J., 2006, MNRAS, 370, 289 

Begelman M. C., Rossi E. M., Armitage P. J., 2008, MNRAS, 387, 1649 

Blandford R. D., Begelman M. C., 1999, MNRAS, 303, LI 

Bondi H., 1952, MNRAS, 112, 195 

Bromm V., Loeb A., 2003, ApJ, 596, 34 

Choi J.-H., Shlosman I., Begelman M. C., 2013, ApJ, 774, 149 

Choi J.-H., Shlosman I., Begelman M. C., 2015, MNRAS, 450, 4411 

Costa T., Sijacki D., Trend M., Haehnelt M. G., 2014, MNRAS, 439, 2146 

Coughlin E. R., Begelman M. C., 2014, ApJ, 781, 82 

Davies M. B., Miller M. C., Bellovary J. M., 2011, ApJ, 740, L42 

Devecchi B., Volonteri M., 2009, ApJ, 694, 302 

Di Matteo T., Colberg J., Springel V., Hernquist L., Sijacki D., 2008, ApJ, 
676, 33 

Dubois Y., Pichon C., Devriendt J., Silk J., Haehnelt M., Kimm T., Slyz A., 
2013, MNRAS, 428, 2885 

Dubois Y, Volonteri M., Silk J., Devriendt J., Slyz A., Teyssier R., 2015, 
preprint, (arXiv: 1504.00818) 

Fan X., et al., 2006, AJ, 131, 1203 

Ferrara A., Salvador! S., Yue B., Schleicher D., 2014, MNRAS, 443, 2410 
Haiman Z., 2013, in Wiklind T., Mobasher B., Bromm V., eds. Astrophysics 
and Space Science Library Vol. 396, Astrophysics and Space Science 
Library, p. 293 (arXiv: 1283.6075), doi:10.1007/978-3-642-32362- 
1_6 

Jeon M., Pawlik A. H., Greif T. H., Glover S. C. O., Bromm V., Milosavl¬ 
jevic M., Klessen R. S., 2012, ApJ, 754, 34 
Jiang Y.-F, Stone J. M., Davis S. W., 2014, ApJ, 796, 106 
Johnson J. L., Bromm V., 2007, MNRAS, 374, 1557 


MNRAS 000, 1-12 (0000) 





12 F. Pacucci et al. 


Johnson J. L., Whalen D. J., Fryer C. L., Li H., 2012, ApJ, 750, 66 

King A., 2003, ApJ, 596, L27 

King A. R., 2010, MNRAS, 402, 1516 

Lasota J.-R, 2015, preprint, (arXiv: 1505.02172) 

Latif M. A., Schleicher D. R. G., Schmidt W., Niemeyer J., 2013a, MNRAS, 
433, 1607 

Latif M. A., Schleicher D. R. G., Schmidt W., Niemeyer J. C., 2013b, MN¬ 
RAS, 436, 2989 

Latif M. A., Niemeyer J. C., Schleicher D. R. G., 2014, MNRAS, 440, 2969 
Lodato G., Natarajan P., 2006, MNRAS, 371, 1813 
Lodato G., Pringle J. E., 2007, MNRAS, 381, 1287 
Madau R, Rees M. J., 2001, ApJ, 551, L27 
Madau R, Haardt F, Dotti M., 2014, ApJ, 784, L38 
McKinney J. C., Tchekhovskoy A., Sadowski A., Narayan R., 2014, MN¬ 
RAS, 441, 3177 

Milosavljevic M., Bromm V., Couch S. M., Oh S. R, 2009, ApJ, 698, 766 
Mineshige S., Kawaguchi T., Takeuchi M., Hayashida K., 2000, PASJ, 52, 
499 

Mortlock D. J., et al., 2011, Nature, 474, 616 

Novak G. S., Ostriker J. R, Ciotti L., 2011, ApJ, 737, 26 

Novak G. S., Ostriker J. R, Ciotti L., 2012, MNRAS, 427, 2734 

Ohsuga K., Mineshige S., Mori M., Umemura M., 2002, ApJ, 574, 315 

Pacucci F., Fen-ara A., 2015, MNRAS, 448, 104 

Pacucci F, Ferrara A., Volonteri M., Dubus G., 2015, preprint, 
(arXiv:1586.05299) 

Paczynski B., Abramowicz M. A., 1982, ApJ, 253, 897 
Park K., Ricotti M., 2012, ApJ, 747, 9 
Petri A., Ferrara A., SalvateiTa R., 2012, MNRAS, 422, 1690 
Planck Collaboration et al., 2015, ArXiv e-prints 1502.01589, 

Sadowski A., 2009, ApJS, 183, 171 
Sadowski A., 2011, ArXiv e-prints 1108.0396, 

Sadowski A., Narayan R., McKinney J. C., Tchekhovskoy A., 2014, MN¬ 
RAS, 439, 503 

Shakura N. L, Sunyaev R. A., 1973, A&A, 24, 337 

Shang C., Bryan G. L., Haiman Z., 2010, MNRAS, 402, 1249 

Silk J., Rees M. L, 1998, A&A, 331, LI 

Springel V., et al., 2005, Nature, 435, 629 

Thorne K. S., 1974, ApJ, 191, 507 

Volonteri M., 2010, Astronomy and Astrophysics Review, 18, 279 

Volonteri M., Rees M. J., 2005, ApJ, 633, 624 

Volonteri M., Silk J., Dubus G., 2014, ArXiv e-prints 1401.3513, 

Wu X.-B., et al., 2015, Nature, 518, 512 


MNRAS 000, 1-12 (0000) 


